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This note is a preliminary account of research undertaken jointly with G. Marmo of Napoli and P. 
Urbanski of Warsaw. 

We propose a new description of dynamics of autonomous mechanical systems which includes the 

C*^ ' momentum-velocity relation. This description is formulated as a variational principle of virtual action 

^^ , more complete than the Hamilton Principle. The inclusion of constraints in this description is the 

main topic of the present note. We give examples and models of constraints in variational formulations 

^ , of statics and dynamics of autonomous systems. 

f^ ' A complete description of the dynamics of a mechanical system involves both external forces and 

C^^ , momenta. In a fixed time interval the dynamics is a relation between the motion of the system in 

^-P ' configuration space, external forces applied to the system during the time interval, and the initial 

and final momenta. This relation is derived from a variational principle which involves variations of 
iQ ' the end points. Constrained systems are idealized representations of unconstrained systems. Such 

f^ , idealizations are appropriate when forces at our disposal are unable to move the configuration of the 

system away from a subset of the configuration space by a perceptible distance. This description fits 
Q.. at least the case of holonomic constraints. We believe that constraints should be imposed on virtual 

displacements. Holonomic constraints are usually interpreted as restrictions on configurations of a 
mechanical system. Nonholonomic constraints are additional restrictions imposed on velocities. This 
C^ ' traditional terminology is not adapted to our concept of constraints as imposed on virtual displacements 

and only indirectly affecting configurations and velocities. Our concept of nonholonomic constraints 
makes perfect sense for static systems even if velocities do not appear in the description of such 
systems. We will use the terms configuration constraints and velocity constraints instead of holonomic 
and nonholonomic constraints. 



2. Geometric structures. 

Let Q be the Euclidean affine space of Newtonian mechanics. The model space for Q is a vector 
space V of dimension 3. The Euclidean structure is represented by a metric tensor g:V ^ V* . The 
space V* is the dual of the model space. The canonical pairing is a bilinear mapping 

(,):l/*xy^R. (1) 

We denote by gi — q^ the vector associated with the points qq and qi. We write qi = qo + v if 
V = qi — qg. The norm ||w|| of a vector u g T^ is defined by 



\\v\\ = ^{gi^vj:^. (2) 

The derivative of a function F:Q^M. is the mapping 

DF:QxV^R:iq,v)^-^Fiq + sv)\s^o- (3) 

ds 
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The first and second derivatives of a difFerentiable curve 7: M ^ Q are mappings 7: R ^ V^ and 
-f.R^V. 

The tangent bundle TQ is identified with Q x V, the cotangent bundle T*Q is identified with 
Q X V* , the second tangent bundle T^Q is identified with Q x V x V, the iterated tangent bundle 
TTQ is identified with Q x V x V x V, and the tangent of the cotangent bundle TT*Q is identified 
with Q X V* X V xV* . We have the projections 

TQ:JQ^Q:{q,q)^q, (4) 

Ttq: TTQ -^ TQ: {q, q, Sq, Sq) 1-^ {q, q), (5) 

Ttq: TTQ -^ TQ: (q, q, Sq, Sq) h^ {q, Sq), (6) 

and the canonical involution 

kq: TTQ -^ TTQ: {q, q, Sq, Sq) ^ {q, Sq, q, Sq). (7) 

For each subset C of Q we have the tangent set 

TC ~ {{q, Sq) G TQ ; there is a curve 7: K ^ Q 

such that 7(0) = g,D7(0) = Sq, and 7(3) e C if s > 0} (8) 

The space 

l-Q = QxV = {{q,v)eJQ;vj^O} (9) 

is the tangent bundle with the zero section removed. 

3. Statics of a material point. 

We consider the statics of a material point in the Euclidean affine space Q of Newtonian physics. 
An element {q, Sq) of TQ is a virtual displacement and an element [q, f) of T*Q represents an external 
force. The evaluation 

{{q,f)Aq,Sq))^{f,Sq) (10) 

of an external force {q, /) € T Q on a virtual displacement {q, Sq) G TQ is the virtual work performed 
by an external device controlling the configuration of the system. 

Admissible displacements form a subset C^ C TQ. If {q,Sq) is an admissible displacement, then 
{q, kSq) is again an admissible displacement for each number fc > 0. The set C^ represents constraints 
imposed on virtual displacements. Implicitly it restricts admissible configurations to the set 

C" = {<? e Q; (q, Sq) e C^ for some SqcV] . (11) 

The inclusion C^ C TC*^ is usually satisfied. We say that constraints are configuration constraints if 
C^ = TC°. The set C — C^ itself is called a configuration constraint. A simple two-sided configuration 
constraint is an embedded submanifold C <Z Q 
There is a function 

cr:Ci->R (12) 

assigning to each admissible virtual displacement the virtual work that an external device has to 

o 

perform in order to effect this displacement. This virtual work function is differentiable on C^ n TQ 
and positive homogeneous in the sense that 

cr{q,kSq) — ka{q,Sq) (13) 

if fc > 0. A typical example of a virtual work function is the mapping 

a:C^ ^R:{q,Sq)^DU{q,Sq) (14) 
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derived from an internal energy function U defined in a domain large enough to make the derivative 
DU{q, Sq) meaningfuL In the case of a configuration constraint C C Q it is enough to have the internal 
energy defined on C. The function 

<j:C'^R:{q,5q)^p{q)\\Sq\\ (15) 

represents virtual work due to friction. 

The response of the system to external control is represented by a set S C T*Q of external forces 
satisfying the principle of virtual work 

{.f,(^<l) < <^{'l,Sq) for each virtual displacement {q,Sq) E C^ . (16) 

The set S is the constitutive set of the system. It can be viewed as the list of possible configurations of 
the system together with external forces compatible with these configurations. If a{q, ~Sq) — —(y{q, Sq) 
and the constraints are two-sided, then the principle of virtual work assumes the simpler form 

(./j Sq) — cr{q, 6q) for each virtual displacement {q, Sq) E C^ . (17) 

Example 1. Let a material point be constrained to a circular hoop with the center a.t qo E Q and 
radius a in the plane orthogonal to a unit vector n E V . We have two-sided configuration constraints 

C° = {geO; {g{q - qo),n) =0,\\q - qo\\ = a} , (18) 

Cl == TCO - {{q, Sq) E TQ; {g{q - qo),n) = 0,\\q-qa\\= a, {g{Sq), n) = 0, {g{Sq), u{q)) = 0} , (19) 
where u{q) is the unit vector {q — qo)\\q — qo\\^^- The constitutive set 

S={{q, /) E T*Q; qEC", (/, Sq) = for each {q, Sq) E C^} 

= {(g,/)GT*Q; {g{q-qo),n)^0,\\q-qo\\=a,{.f,u{q))^0,{f,n)^0} (20) 

represents the statics of the system without friction and the constitutive set 

S^{{q,f)ET*Q;qEC°,{f,Sq)<p\\Sq\\ for each {q,Sq)EC^} 

- {{q, I) e T*g; {g{q - qo),n) = 0, \\q - qo\\ = a, {f,u{q))^ + (/, nf < p^} (21) 

takes constant friction into account. ▲ 

Example 2. Let a material point be constrained to the exterior of a solid ball with the centre at 
qo E Q and radius a. In its displacements on the surface of the ball the point encounters friction 
proportional to the component of the external force pressing the point against the surface. Correct 
representation of the statics of the point is obtained with one-sided constraints 

C'' = {qEQ; \\q-qo\\>a}, (22) 

C = {(g, Sq) E TQ; \\q - qo\\ > a, {g{Sq), u{q)) > i^^\\Sq\\^ - {g{Sq), u{q))^ if \\q - qo\\ = a] (23) 
and the constitutive set 

S^{{q, f) E T*Q; q E C°, (/, Sq) < for each (g, Sq) E C^} 
= {(«,/) eT*Q; ||(7-go||>a,/ = if \\q - qo\\ > a, 

i^{f,u{q)) + ^\\f\\^-{f,u{q))^<0 if llrj-goll-a}, (24) 

where u{q) = {q — Qo)\\<l ^ ^oll"^- The constraints in this example are not configuration constraints. ▲ 



Example 3. Let i, j, and k be mutually orthogonal unit vectors and let qo be a point. Let one-sided 
configuration constraints be specified by 



C" = {qe Q; {g{q -qo),l)>Q, {giq - qo),j) > 0} (25) 



and 



C^ = TC° 

= {{q, Sq) e TQ; {g{q - qo), i) > 0, {g{q - go), j) > 0, 

{g{dq),i)>0, if(5(g-go),*) = 0,(g(J<z),j}>0, if (g(g - <Zo), j) - 0} . (26) 

The statics of a material point not subject to internal forces is represented by the constitutive set 

S^{{q,f)eT*Q;qeC",{f,Sq)< for each {q,Sq)eC'} 
= { (9, /) e T*Q; {giq - qo), i) > 0, {g{q - qo),j) > 0, (/, k) = 0, 

(/,z)=0 and {f,j)<0 if {g{q-qo),j)^0 and {g{q - qo),i) ^ 0, 
{f,j)=0 and (/,*)< if {g{q-qo),i)=0 and {g{q - qo),j) ^ 0, 
{f,i) < and (/,j) < if {giq - qo),z) = and {giq - qo),j) = 0} (27) 



Example 4. In terms of the vectors i, j, and k and the point qo of the preceding example we define 
one-sided configuration constraints by 

C° = {qeQ; {giq-qo),l)<0 or {giq - qo),j) < 0} (28) 

and 

C^ = TC" 

= {iq,Sq)eTQ; {giq - qo),i) < or (g(g - go),j) < 0, 

{gi<>q),i)<0 if (.g(g-(Zo),») =0 and (.9((7 - go), j) ^ 0, 
(5('5'7),j) < if (.g('?-'Zo),i) ==0 and (.g(g - go),«) # 0, 
{giSq),i)<0 or (g(<5g),i)<0 if (^(g - go),*) = and (g(g - go), j) = 0} . (29) 

The statics of a material point not subject to internal forces is represented by the constitutive set 

5-{(g,/)eT*Q; geC",(/,5g)< for each iq,6q)eC'} 
= {(g,/)eT*Q; (,g(g-go),«)<0 or (.g(g - go),.?) < 0, (/, fc) = 

{.f,t)=0 and (/,j)>0 if (.g(g-go),.?)-0 and {giq - qo),t) ^ 0, 
(/,.?) =0 and (/,*)> if (.g(g-go),*)-0 and (^(g - go),.?) # 0, 
(/,^) > and (/,j) > if {giq - qo),i) = and {giq - go),j) - 0} (30) 



4. Modeling configuration constraints in statics. 

We believe that constraint static systems are idealized representations of unconstrained systems. The 
magnitude of the force that an external device can apply to a static system is limited and instruments 
used to observe displacements have a limited resolution. Idealizations take these limitations into 
account. We restrict the analysis to configuration constraints. A definition of configuration constraints 
will be based on the assumption that the norms of external forces at our disposal have an upper bound 
F and that displacements of distances less than d can not be detected. 



Let C be a subset of Q. We denote by d{q, C) the distance of a point q E Q from C. If C C Q 
is an embedded submanifold or a submanifold with smooth boundary, then for each configuration q 
sufficiently close to C there is a unique point qc E C nearest to q. The distance d{q,C) ~ \\q — qc\\ 
of q from C is a well defined function in a neighbourhood of C. If q is not in C, then qc ^ q and the 
unit vector e{q) = {q — qc)\\<l ^ QcW^^ is orthogonal to C or the boundary of C at qc- 

Let C C Q be an embedded submanifold or a submanifold with smooth boundary and let 5* C T*(5 
be the constitutive set of a static system derived from the principle of virtual work 

(/, 5q) < (7{q, Sq) for each virtual displacement {q, 6q) G TC. (31) 

A model of this static system is constructed by choosing a function a on TQ such that a is the 
restriction of a to TC and replacing the original principle of virtual work by the principle 

(/, Sq) < a(q, 6q) + kd{q, C){g{e{q)), Sq) for each virtual displacement {q, Sq). (32) 

The term kd{q,C){g{e{q)),6q) is the directional derivative DK{q,Sq) of the elastic internal energy 
function 

Kiq)^'^id{q,C))' (33) 

defined in the neighbourhood of C in which the distance function d{q, C) is well defined. The inequal- 
ities 

kd{q,C)<{f,e{q))+a{q,-e{q)) (34) 

and 

kd{q,C)<\{f,e{q))\ + \a{q,-e{q))\ (35) 

are derived from the principle of virtual work by setting Sq = —e{q). We will assume that F <S^ kd and 
expect that the inequahty \<T{q, —e{q))\ <C kd is satisfied. These inequalities together with \{f,e{q))\ < 
11/11 < F result in d{q, C) ^ d. It follows that using external forces at our disposal we can not induce 
the material point to assume configurations at noticeable distances away from C. It also follows that 
within the limits imposed by ||/|| < F the component {f,e{q)) is arbitrary. Examples will be used to 
clarify details and present variations of this construction. 

Example 5. Let C be the set C° of Example 1. We obtain the equation 

{d{q, C))2 = \\q - qof - 2aVlk - 9oP - {giq - qo),n)^ + a' (36) 

for the distance d{q, C) of q from C if this distance is less than a. If (i(g, C) ^ 0, then 

e{q) = g-go-(.9(g-go)^")" (37) 

Vlk-9o|p - {g(q-qa),n)^ 

is the unit vector orthogonal to C at the point q' E C closest to q pointing from q' to q. Let a function 
a:TQ ^ M be defined by a{q,n) — 0, a{q,q — qo) = 0, and a{q,Sq) = p\\Sq\\ if {g{Sq),n) = and 
{g{Sq), q — qo) = 0. The unconstrained system represented by the principle of virtual work 

(/, Sq) < CT(g, Sq) + kd{q, C){g{e{q)), Sq) for each virtual displacement {q, Sq) (38) 

is a model of the constrained system of Example 1. It follows from the principle of virtual work that 
/ = kd{q,C)g{e{q)) + /', where the component /' satisfies relations {f',n) = 0, {f',q — qo) — 0, and 
ll/'ll ^ P- If fc ^ oo, then d{q, C) -^ 0. Any value can be obtained for the component kd{q, C)g{e{q)) 
as fc ^ oo and d{q, C) -^ 0. This is in agreement with the principle of virtual work of Example 1. ▲ 



Example 6. Let C be the set C° of Example 2 and let u{q) = (q — qa)\\q — qo\\ ^- The distance 
d{q, C) is equal to \\q — q^yW — a. A function a: TQ ^ M is defined by (T(g, 5q) = if \\q — <Zo|| > a and 

a{q, 6q) ^ -kd{q, C){g{u{q)),Sq) + kd{q, C)v^\\5qr-{g{u{q)),5q)^ (39) 

if \\q — <7o|l < a. The principle of virtual work 

(/i 5q) < <^{<ii Sq) for each virtual displacement (g, 5q) G TQ (40) 

implies the following relations for the external force /. If \\q — go II > «: then / = 0. If ||g — QoH < a, 
then / = -kd(q,C)g{u{q)) + r with {f ,u{q)) = and (/',,5(z) < kd{q,C)u\\5q\\ if ((7(m((z)), <5g) - 0. 
If fc ^ oo, then d{q,C) -^ 0. The component {f,u{q)) ~ —kd(q,C) can have any negative limit and 
(/', 5q) < — (/, u{q))iy\\Sq\\. This is in agreement with the principle of virtual work of Example 2. ▲ 

Example 7. A model for the system in Example 3 can be easily constructed even if the boundary 
of the set C = C" is not smooth. The distance d{q, C) is defined by 

diq,C) = -{g{t),q-qo) (41) 

if {9{j),Q-Qa) > and {g{i),q-qo) < 0, 

d{q,C)^-{giJ),q~qo) (42) 

if (5(0; 9 - Qo) > and {g{j),q- qo) < 0, and 



d{q, C) - v/(.9(*),'7-'Zo)2 + (.g(.7),g-go)2 (43) 

if {g{i),q — qo) < and {g{j),q — qo) < 0. A vector field e{q) is defined is defined outside of C by 
e{q) = -i if {g{j), q - qo) > and {g{i),q - qo) < 0, 

e{q) = -J (44) 

if {g{i),q - qo) >0 and (.g(j), q - qo) < 0, and 

e(g) = Ml), q - qo)l + {gij),q - qo)j){diq, C))-' (45) 

if {g{i), q — qo) < and {g{j), q — qo) < 0. A function a on TQ is defined by a{q, dq) = if g G C and 
(7{q, dq) = k{g{e{q)), 6q) \i q ^ C . The constitutive set of Example 3 is obtained from the principle of 
virtual work 

ifi^Q) — '^{Qt^q) foJ' each virtual displacement {q,5q) G TQ (46) 

with k ~> oo and d{q, C) ^ 0. ▲ 

Example 8. The construction of the model in the preceding example followed exactly the prescription 
given at the beginning of the present section. This construction can not be directly applied to the set 
C = C'^ of Example 4. It can be applied to the modified set 

Cs = Co\ {q G Q; {g(i),q - qo) < r, {g{j),q - qo) < r, 

{{9{^),q - qo) - r)^ + ((^(j), 9 - ^o) - r)^ > r} . (47) 

The original set Co is obtained as the limit as r ^ 0. ▲ 

6 



5. Kinematics of autonomous systems and scleronomic constraints. 

Motions of a material point in the Euclidean afRne space Q are curves ^: I ^ Q parameterized by 
time t in an open interval / C M. We have the tangent prolongation (^, ^): / ^ TQ: t h^ {(.(t), ^{t)) and 
the second tangent prolongation {£.,£,,£,)■ I -^ T^Q:t i-^ ii(t),£,{t),£,{t)) of a motion ^. 

Variational formulations of analytical mechanics require the concept of a virtual displacement of a 
motion. A virtual displacement of a motion ^ is a mapping {^,SS,)'.I -^ TQ:t h^ (^(t), (5^(t)). This 
mapping is obtained from a homotopy 

X:RxI^Q. (48) 

The base curve x(0, •) is the motion ^. The virtual displacement is the mapping 

(e,<5e):/^TQ:t^tx(-,t)(0). (49) 

A mapping (^, ^, S^, S^): I — > TTQ: t h^ {(.{t),^{t),6^{t), S(,{t)) is obtained from a virtual displacement 
(^, 6^) as the composition Hq o (^, S^, ^, 6^) of the tangent prolongation (^, 6^, ^, S^) with the involution 
kq. Virtual displacements are subject to constraints. All considered versions of constraints can 
eventually be reduced to differential equations formulated in terms of a subset C'^'^^ C TTQ such that 
if {q,q,Sq,Sq) £ C'-^'^' , then {q,q,k5q,kSq) G C'-^'^' for each number A; > 0. An admissible virtual 
displacement (^, 6^) is required to satisfy the condition 

{m.m.sm^smec-'^'^ (50) 

for each t ^ I . This condition implies conditions 

{m.m)ec^'''\ (51) 

(C(i),<5e(t))eC(i-"), (52) 

and 

e(t) e C(°'") (53) 

for each t e /. Sets C^^'^\ C^^'^\ and C("^°) are defined by 

C(°'l) = {{qA) e Tg; {q,q,5q,5q) E C^^'^) for some {dq,6q) e V X v} , (54) 

C(l'°) = I {q, Sq) e jg; (g, q, Sq, Sq) E C^^^^^ for some {q,6q)EVxV^, (55) 

and 

C("^°) == L e g; (g, ^) e C(°'i) for some q E v\ . (56) 

Condition (51) is a differential equation for the motion ^: I ^ Q. Constraints arc usually discussed 
it terms of this equation. The inclusion 

^(0,1) ^ jc(o,o) (57) 

must be satisfied since it is a necessary integrability condition for the equation (51). Constraints will 
be called configuration constraints if 

C(o,i)^TC<ofi), (58) 

Velocity constraints are said to be ^mear if the set C^^''^^ is a submanifold of g and C^^'^^ is a distribution 
on this submanifold. Linear constraints are said to be holonomic if C^^'^' is integrable in the sense 
of Frobenius. Configuration constraints are a special case of holonomic constraints. Sets C^^'^' and 
C^^'^' are not usually discussed directly even if information contained in the velocity constraints C^'^'^' 
is not sufficient for the application of variational methods. The condition (50) is equivalent to 

im,mt), m, sm) - t(e, soit) e «q(c(1'1)) c jjg. (59) 
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It is a differential equation for tlie virtual displacement (^, S^): I -^ TQ. The inclusion 

C(i'i) c kq(TC(i^o)) (60) 

is a necessary integrability condition for this equation. 

Two different methods of constructing the set C^^'^^ from the velocity constraints C^^'^^ are found 
in an article of Arnold, Kozlov, and Neishtadt [Arn] . 

(1) In vaconomic mechanics the natural construction 

^(1,1)^7^(0,1) (61) 

is used. This construction is the result of the differential equation 

tx(s,.)(t)eC(°'i) (62) 

imposed on curves x(s, •) also with s ^ 0. See [Arn] for modifications of this construction 
necessary when virtual displacements vanishing at the ends of a time interval are required. 
With these modifications the formula (61) is still valid. The set C'^^'"-' is the tangent set TC^'^'^^ 
ofC(0'"). 

(2) The d'Alembert- Lagrange principle is based on the inclusion 

C(l'l) C C(mT ^ I (g, 5, Sq, dq) G TTQ; {q, q, 0, dq) G TC^"^'^ } . (63) 

This construction derives from the condition 

tx(-i)WeC(i'0) (64) 

for each t G / and each s and the condition 

tx(0,.)(t)eC("'i) (65) 



imposed on the base curve S, — x(Oi ') but not the curves x(s, •) for s ^ 0. The set C^^'^^ may 
not represent an integrable differential equation (59) for (^, 5S,)- The set C^^'^^ is the integrable 
part of C(i'i). If C(°'i) is a vector subbundle of TQ, then C^^'"^ = C^"-^\ If C^°''^^ is an affine 
subbundle, then (7*^^'°' is the model bundle. In both cases 

C(My ^ I (g, q, 6q, 6q) G TJQ; {q, q) G C^''''^ , {q, Sq) G ^(1'°) } (66) 

and 



^(1,1) ^ c(^A) n T(^^^(c(i,i))), (67) 

For configuration constraints both construction give the same result 

C(i'i)=TTC(°'"). (68) 

6. Dynamics of unconstrained autonomous systems. 

We see four possible formulations of dynamics. 

A. Dynamics of a material point can be specified as a collection of boundary value relations with 
external forces associated with time intervals. A boundary value relation for a time interval [a, 6] C M 
is a relation between an arc ^: [a, b] ^ Q, a, mapping (^, ip): [a, b] -^ T*Q, and two covectors (C(a), 7r(a)) 
and (^(6), 7r(&)). The mapping {S,,ip) represents the external force applied to the material point along 
the arc ^ the covectors (^(a),7r(a)), (^(fe),7r(6)) are the initial momentum and final momentum. It 



is convenient to consider the arc ^: [a,b] -^ Q and the mapping (C^p)- [a,b] -^ T*Q the restrictions 
to the interval [a,b] of mappings (,'■ I ^ Q and (^,(^):/ -^ T*Q defined on an open interval / C M 
containing [a,b]. The covectors (^(a),7r(a)) and (^(&),7r(6)) will be considered values of a mapping 
{£,, it): I -^ T*Q at the ends of the interval. The two mappings {£,, (fi): I -^ T*Q and {£,, n): I -^ T*Q 
can be combined in a single mapping 

{t^,TTy.I^QxV* xV*. (69) 

An element 

((e, ^): [a, b] ^ T*Q, (^(a), 7r(a)), (^(6), 7r(6))) (70) 

of the boundary value relation Dt^u for an interval [a, b] satisfies the virtual action principle 

' {^{t),sat))dt - {nib), sm) + (^(«), m^)) 

' Km,m,mt),sm) - m(gm),sm) dt. (71) 

for each virtual displacement (^, 5S,)'. [a, b] -^ TQ obtained as a restriction to [a, b] of a virtual displace- 
ment (C, 6^): I -^ TQ. The term 

m{gm),Sm (72) 

is the derivative DT{^{t),^{t),S^{t),6^{t)) of the kinetic energy function 

T:TQ^R:{q,q)^j\\qf (73) 

of a material point with mass m. The function A: TTQ -^ K represents the virtual action of internal 
forces. For the sake of simplicity we assume that A is a linear form (a linear function of {5q, Sq)). The 
proposed principle of virtual action is more general than the Hamilton Principle. Note that the virtual 
displacements {^{a),d^{a)) and {£,{b) , S£,{b)) of the end points of the arc do not vanish. Variational 
principles with variations of end points but without external forces were considered by Schwinger. The 
momentum- velocity relation is a law of physics and is a part of dynamics of a material point. This 
relation is included in the Schwinger version of the principle of virtual action but not in the Hamilton 
Principle. 

Examples of the virtual action function A include the function 

A: TTQ ^ M: (q, q, Sq, 5q) ^ e{A{q), Sq) + e{DA{q, Sq),q) (74) 

for a charged particle in a magnetic field derived from the vector potential A:Q ^ V* and the function 

A: TTQ ^ M: (g, q, Sq, 5q) ^ j{g{q), Sq) (75) 

for a material point immersed in a viscous medium. In the first of these examples the function 
X{q, q, Sq, Sq) is the derivative Da{q, q, Sq, Sq) of the function a: TQ -^ R: {q, q) ^^ e{A{q),q). We will 
continue the analysis assuming that the function A is of the simpler type 

A: TTQ ^ M: (q, q, Sq, Sq) ^ {fi{q),Sq), (76) 

where /i is a mapping from TQ to V* . The principle of virtual action assumes the simpler form 

b 

{^{t), Sm)dt - {n{b),Sm) + ^(a), 5e(a)) 



{{i^im),sm) - migim, sm) a. (77) 
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Equivalent versions of this variational principle 

rb 



{ip{t)Mt))dt - (7r(6), Sm) + (^(«), Saa)) 



m{g{i{b)), Sm) + miaiiia)), S^a)) (78) 



and 



b 



{wm\ sm) - Mgim, sm) a (79) 



are easily derived by using the identities 

' 7r(t), sm) + {At), sm)) dt= I ±{n{t),6m)dt 



{7r{b),dm))-{Aa),Sm)) (80) 



and 



m{gim), sm)) + Maim, sm )dt= 1 mj^igim, sm)dt 



^ m{g{m), dm)) - m{9im)), sm))- (81) 



B. Dynamics can be specified as the collection D of curves 

{^,^,n):I^QxV* xV* (82) 

defined on open intervals / C M with the property that for each time interval [a, b] d I the arc 
(^, </?)|[a, 6] and the covcctors (^(a),7r(a)) and (^(6),7r(6)) arc in the boundary relation -D[a,6]- 



C. Dynamics can be specified as differential equations 

^(t) - m) = f^im) (83) 

and 

m = giiit)) (84) 

for mappings (^, (/3, tt): / ^ Q x V* x V* . These equations will be denoted by D. 



D. Dynamics can be specified as differential equations 

^{t) ^ mgim) + Km) (85) 

and 

<t) ^ gim (86) 
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for mappings {^,(p,tt):I -^ Q x V* x V* . These equations will be denoted by E. 



Of the four formulations of dynamics version A is fundamental. The family of curves 2) and the 
differential equations D and E introduced in B, C, and D are auxiliary objects. 

Differential equations D and E are obviously equivalent. 

We show that the family CD is the set of solutions of the equations D. The proof is based on version 
(79) of the principle of virtual action. If a curve (^, (/3, tt): / ^ Q x V* x V* is in D and {£,, SS,): I -^ TQ 
is an arbitrary virtual displacement, then the equality (79) holds for all intervals [a, b] C /. It follows 
that the equality 

Mt) - Mt), sat)) - {7r{t), sm - {^l{at)),sm) - migim, sm (87) 

holds at each t E I. This implies that equations D are satisfied due to arbitrariness of the vectors S^t) 
and S^t)- Conversely if (^, (^,7r):/ ^ Q x V* x V* is a solution of D, then the equality (87) holds 
in I with an arbitrary displacement (^, 6^): I —>■ TQ. The validity of the principle of virtual action for 
each time interval [a, 6] C / is established by integration. Hence, (^, (p, n) is in 2). 

It follows from the definition of CD that this family is constructed from the boundary value relations. 
We show that elements of boundary value relations can be constructed from elements of CD. Let [a, b] 
be a time interval included in an open interval / C K. and let (^, (/3, tt): / ^ Q x V* x V* be a mapping 
such that (^, ip)\[a, b], (^(a), 7r(a)), (^(&), 7r(6))) is in £'[a,6]- It follows from version (78) of the principle 
of virtual action that the mapping (^, (p) satisfies the equation 

fit) = rng(m) + mt)) (88) 

in [a, b] and that 

7:{a) = g{i{a)) and 7:{b) - g(^(5)). (89) 

Let mappings ip': I -^ V* and tt': / -^ V* be defined by 

f'it)^mgm) + ^^im) (90) 

and 

7:'{t) - mgiiit)). (91) 

The mapping (^, (/3',7r'): / ^ Q x V* x V* is in D since it is a solution of E. The boundary value 
data extracted from this mapping are in the boundary value relation since (^, (/?')![&, 6] — (^, (y9)|[a, 6], 
7r'(a) = 7r(a), and 7r'(6) = 7r(6). 

7. Dynamics with configuration constraints. 

We list four possible formulations of dynamics with constraints analogous to the four formulations 
in the preceding section. As defined in Section 5 an admissible virtual displacement is a mapping 
(^, (5^): / -^ TQ satisfying the condition 

{m,i{t),Sm, Siit)) e C(i'i) = TTC(°'") (92) 

for each t E I. 

A. Dynamics of a material point can be considered a collection of boundary value relations associated 
with time intervals. An element 

((e, f): [a, b] ^ 7*0, (e(a), 7r(a)), (^(6), n{b))) (93) 

of the boundary value relation D^a.b] for an interval [a, b] satisfies the virtual action principle 
b 
{^{t),Sat))dt - {n{b),Sab)) + (7r(a), <5e(«)) 



{wm), sm) - Mgim, sm) dt. m 
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for each admissible virtual displacement (^,5^): [a,b] -^ TQ obtained as a restriction to [a, 6] of an 
admissible virtual displacement {S,,5£_):I -^ TQ. The mapping /i is defined on C^^'"^ = TC'°'"'. The 
condition ^(i) £ (7*^"'°^ for each t e / is implied. 

There are again the equivalent versions of this variational principle 

b 

{ip{t),sat))dt - {nib), sm) + (^(«), ma)) 

' m{9m),sm) + {i^im),sm)] dt 



and 



migiiib)), Sm) + miaUia)), S^a)) (95) 



' {^(t) - 7:it),Sm) - {At),Sm) dt 



{{i^im),sm) - Mgim, sm) dt. (96) 



B. Dynamics can be specified as the collection D of curves 

{^,ip,7r):I^QxV*xV* (97) 

defined on open intervals / C M such that for each time interval [a, b] C I the arc (^, (f)\[a, b] and the 
covectors (^(a),7r(a)) and (^(6),7r(&)) are in the boundary relation D^^^. 



C. Dynamics can be specified as the differential equation 

(^(i) - n{t), sm) {n{t), 5i{t)) = (^(e(t)), sm) - {gim, sm (98) 

to be satisfied by a curve {^,(p,tt): I ^ Q x V* x V* at each t E I and for each (q, ^(t), S^{t), S(,{t)) G 
(7(1.1) _ ^j^jg jg equivalent to equations 

{ip{t)-7r{t),m))^{mt)),sm) m 

and 

{m,m))^{9m),sm) (loo) 

satisfied at each t e / for each (^(t), S^) ^ C^'^'°'>. 



D. Dynamics can be specified as the differential equations 

{^it),5m) - {mgm) + ^iim), sm) (loi) 

and 

{m),sm)^{9m),sm) (102) 

satisfied by a mapping (^, lp,tt): I ^ Q x V* x V* at each t e / and each (^(t), Sm) ^ C'-^^'^K 



The four formulations are valid for configuration constraints and are equivalent as in the case of 
unconstrained systems. The situation is much more complex in the case of more general constraints. 
The method of models could be a tool for testing the validity of different formulations. We will apply 
this tool to the momentum-velocity relation. Note that the usual momentum- velocity relation 

m) = gim (103) 

is replaced by the equation (102). We will attempt a justification of this modification of the momentum- 
velocity relation based on models of configuration constraints. 
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8. Models of autonomous systems with configuration constraints. 

Let a material point of mass m be constrained to a plane C^^'^' C Q passing through a point qo 
and orthogonal to a unit vector n. We will assume that there are no internal forces and no external 
forces are applied. The constraint will be modeled by a strong internal elastic force k{g{q — qo), n)g{n). 
Let an initial momentum (^(a),7r(a)) such that ^(a) G C^^'^^ and (7r(a),n) = be applied to the 
point. The solution of the dynamical equations will be the mapping {£^,ip,TT): R ^ Q x V* x V* with 
^(i) = ^(a) + m~^g~^{iT{a)){t — a), (p{t) = 0, and 7r(i) = 7r(a). If the constraint is replaced by the 
elastic force the solution mapping will be the same. Let now the initial momentum have a non zero 
component (7r(a),n). For the unconstrained model the solution is the mapping (^, (/?, tt) with 

C(t) =C(a) +m^^ {g^^{TT{a)) - (7r(a),n}n) (t - a) + w"^(7r(a), n)nsinw(t - a), (104) 

ip{t) = 0, 

7r(t) = n{a) + {n {a), n)g{n) {cos uj{t — a) — 1), (105) 



and ui = y'kjrn. The oscillation may be invisible since the amplitude uj^^{Tr{a),n) may be small due 
to the high value of lo. The rapidly changing component (7r(a),n}g(n) cosuj{t — a) of the momentum 
transverse to the plane C^^'^' depends on the initial value {Tr{a),n)g{n) and is arbitrary within certain 
limits. This component can be detected by making the material point collide with an unconstrained 
mass. Time dependent external forces and curvature of constraint set C^^'^' may even cause the 
transverse component of momentum influence the visible part of the motion along the constraint. 
The element of the boundary value relation for the idealized constrained system is composed of the 
mapping {£_, (p): [a, b] ^ Q x V* with S^(t) = £_{a) + m^^ (^g^^{TT{a)) — (7r(a), n)nj (t — a), if{t) = 0, and 
covectors (^(a),7r(a)) and (^(6),7r(6)) satisfying the equality 7r(6) — (7r(6), n)5(n) — ■K{a) — {'K{a),n)g{n). 
The transverse component (7r(a),n)g(n) of the initial momentum is arbitrary. Due to the rapidity of 
oscillations the final value of the transverse component (7r(6),n)g(n) of final momentum should be 
considered arbitrary and independent of the initial value. 

This analysis based on a purely elastic model suggests that the virtual action principle 

{ip{t)M(t))dt - {Ab),5m) + (^(«),<5C(a)) 



- m{gU{b)),Sab)) + m{gU{a)),Saa)) (106) 

for each admissible virtual displacement (^, S^): [a, b] -^ TQ is appropriate for material points subject 
to configuration constraints. If this formulation of dynamics with configuration constraints is adopted, 
then the momentum-velocity relation 

momentum ~ mass x velocity 
is no longer valid. Velocity is the rate of change of configuration. The component of velocity 
transverse to the constraint set is zero. This is not true of the transverse component of momentum. 
Along directions tangent to the constraint the usual momentum- velocity relation holds. Other models 
may be found appropriate in certain situations. When a rigid bar is struck with a hammer it starts an 
invisible vibration detectable through the sound it emits. The sound is due to momentum (and energy) 
transfer to air molecules colliding with the surface of the bar. The vibration will eventually die away 
although it may continue for a long time if the bar is placed in vacuum. The inevitable damping can be 
taken into account by supplementing the elastic force k(g{q — qo), n)g{n) with a viscous damping force 
7(5(^(0); '^)ff("-)- It may be appropriate to consider the transverse component of the initial momentum 
arbitrary and the transverse component of the final momentum effectively reduced to zero. In cases 
of relatively short time intervals the effect of the damping can be ignored. In cases of strong damping 
and relatively long time intervals. It may be correct to assume that the transverse component of the 
initial momentum is completely absorbed in an essentially inelastic collision with the suspension and 
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not transferred to the material point. In such cases the boundary value relation will be a solution of 
the Hamilton principle with external forces: 



b fb 



{^{t),5m)dt - / (m(5(e(t)), sm) + imt)), sm)) dt (lo?) 

for each admissible virtual displacement {^,SS^):[a,b] -^ TQ with S^{a) = and 6^{b) — 0. The 
equalities 7r(a) = g(^(a)) and 7r(fo) — !j{£,{b)) supplement the variational principle. 
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